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Abstract 


Susceptibilities of conserved charges such as baryon minus lepton number enter baryo- 
genesis computations, since they provide the relationship between conserved charges 
and chemical potentials. Their next-to-leading order corrections are of order g, where 
g is a generic Standard Model coupling. They are due to soft Higgs boson exchange, 


and have been calculated recently, together with some order g 2 corrections. Here we 


compute the complete g 2 contributions. Close to the electroweak crossover the soft 
Higgs contribution is of order g 2 , and is determined by the non-perturbative physics 
at the magnetic screening scale. 
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1 Introduction 


In the early Universe all charges which are violated at a rate smaller than the Hubble 
expansion rate can be considered conserved. For instance, in the minimal Standard 
Model (with zero neutrino masses) baryon number B and the n f = 3 flavor lepton 
numbers Lj are conserved below the electroweak scale, while at higher temperatures 
only the differences X t = B/rif — are conserved. All equilibrium properties are 
determined by the temperature T together with the values of all conserved charges 
Qi or equivalently by the corresponding chemical potentials /ip These properties are 
encoded in the grand canonical partition function 

exp(—fi/T) = trexp [{mQi ~ H)/T ], (1.1) 


where H is the Hamiltonian. 

It is rather plausible that initially the values of conserved charges were practically 
zero, for example if one assumes that the Universe underwent an early period of infla¬ 
tion. Since there is something rather than nothing, some processes must have created 
at least the charge that we know is non-vanishing at present, i.e., the baryon num¬ 
ber, or baryon asymmetry of the Universe. Such a process, called baryogenesis, must 
proceed out of thermal equilibrium. For example, in leptogenesis [T] a non-vanishing 
value of some X, is generated. Afterwards this quantity is conserved and its value de¬ 
termines the equilibrium properties, such as the expectation values of baryon number 
B or lepton number L. 

The values of the charges and thus of the chemical potentials are usually small, so 
that the grand canonical potential is only needed to lowest non-trivial order, which is 
0(/r 2 ) I Then the ft-dependence is fully determined by the second derivatives at zero 
fi, the so-called susceptibilities 


Xi 


i d 2 n 

V dfiidfij 


fi =o 


( 1 . 2 ) 


One important use of the grand canonical potential is to determine the relation 
between B or L and the Qi. Strictly speaking one cannot introduce a chemical potential 
for B + L in the symmetric phase where electroweak sphalerons rapidly violate B + L. 
Nevertheless, one can formally introduce a chemical potential for B + L as long as one 
computes only the expectation value of B + L and not higher moments. The reason is 

3 We assume that the charges Qi are odd under CPT. Then their expectation values vanish when 
y = 0, and fl contains no terms linear in y. 
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that for the resulting partition function 


exp(—Si'/T) = trexp | + L) + fj-iQi - II]iT j 


(1.3) 


one only needs the expansion to first order in [i b +l- Then, even though B + L does 
not commute with H , the operator ordering does not matter because of the trace. The 
expectation value can then be written as 



(1.4) 


This relation can be used to determine B + L and thus B from the value of B — L before 


the electroweak crossover, neglecting possible effects of the non-equilibrium epoch when 
the electroweak sphaleron transitions are shut off. 

Another use of the susceptibilities (1 1.2 [) has been pointed out recently j2] in the 
context of leptogenesis. There the asymmetry can be obtained from a set of kinetic 
equations. One coefficient in these equations quantifies the amount of washout of the 
asymmetry, ft was found that at leading order in the right handed neutrino Yukawa 
couplings the washout rate can be factorized into a product of a spectral function 
which contains dynamical information, and the inverse of a matrix of susceptibilities. 
The spectral function has been computed at next-to-leading order which is 0(g 2 ) in 
the Standard Model couplings It turned out that deep in the symmetric phase the 
NLO corrections to the susceptibilities already start at order g. The 0(g) contribution 
computed in [2j is an infrared effect caused by the exchange of a soft Higgs boson. Close 
to the electroweak crossover the effective thermal Higgs mass can become very small. 
If it becomes of the order of the magnetic screening scale, the perturbative expansion 
for the susceptibilities can be expected to break down. 

In this paper we compute the complete 0(g 2 ) corrections to the susceptibilities, 
thereby completing the 0(g 2 ) result for the washout rate. We obtain contributions both 
from hard (~ T ) and smaller momenta, which, depending on the value of the thermal 
Higgs mass, can be soft (~ gT ) or even smaller (‘ultrasoft’). We use dimensional 
reduction, a framework which allows us to systematically treat the contributions at 
the different scales and the required resummations. 

Part of the 0(g 2 ) susceptibilities have already been computed in |2]. Dimensional 
reduction in the presence of chemical potentials has been considered in |3], where the 

4 For our power counting we make no distinction between the different Standard Model couplings. 
In this respect we differ from [3]. 
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focus was on a electroweak phase transition. Therefore only those terms which depend 
on the Higgs field were computed. 

This paper is organized as follows. In section [2] we recall the role of gauge charges 
and gauge fields in the presence of chemical potentials for global charges. Section [3] 
outlines our use of dimensional reduction. The hard Higgs contribution is obtained in 
section [4] and the dimensionally reduced theory is described in section [5} Depending 
on the value of the effective Higgs mass we obtain either soft (section [6]) or both soft 
and ultrasoft contributions (section [7]). Finally, in section [8] we illustrate our results 
by computing the relation of B and B — L near the electroweak crossover. 


2 Chemical potentials and gauge charges 


We write the partition function (II.ip as a path integral with imaginary time t = —it, 


exp(— VL/T) 


£^<f>exp 



l^iQi + 



( 2 . 1 ) 


where $ stands for all fields in our theory with the Lagrangian «5f. The temporal 
component of the gauge fields act as Lagrange multipliers which enforce Gauss’ law. 
We work in a finite volume and take the volume to infinity in the end. Then, with 
spatial periodic boundary conditions, the total gauge charges vanish. These conditions 
are enforced by the constant modes of the temporal component of the gauge fields. 

In the presence of chemical potentials for global charges the temporal components 
of the gauge fields can develop constant expectation values which act like chemical 
potentials for the corresponding gauge charges. We will only consider the symmetric 
phase of the electroweak theory, where only the weak hypercharge gauge field B M can 
develop an expectation value. 

It is convenient to perform the path integral (12.ip in two steps [0]. First one 
integrates over all fields except over the constant mode of B 0 which we denote by B 0 . 
We denote the result of this integration by exp(— Q/T). In the presence of chemical 
potentials D may contain terms linear in B 0 . The linear terms can arise when some of 
the global charges are correlated with the hypercharge. Then the integral over Bo 


exp(—fl/T) = J dB 0 exp(—n/T) (2.2) 

can lead to //-dependent contributions. 

Here we are interested in small values of the conserved charges which corresponds to 
small values of the chemical potentials. Therefore we need to keep only those terms in 
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O which are at most quadratic in the chemical potentials. Then (12. 2 j) can be evaluated 
in the saddle point approximation, 


exp(— VL/T) = const x exp —12(saddle point )/T 


(2.3) 


Here fl is evaluated at the saddle point 


on 

dB~o 


0, 


(2.4) 


and the constant in (12.31) is independent of the chemical potentials. The relation 
(12.41) determines the expectation value of B 0 and is usually referred to as ‘equilibrium 
condition’. Note that it follows from the saddle point approximation to (12.2ft . 

Our convention is such that the hypercharge gauge field enters the covariant time 
derivative for species a with hypercharge y a as follows, 


D 0 = d t + iy a giB 0 H-= i(d T + y a giB 0 ) H-, (2.5) 

where y v = 1/2 for the Higgs field, and g\ is the weak hypercharge gauge coupling. 
Note that Bq is purely imaginary. The constant mode acts like a chemical potential 
/i a = y a fiY for each species a with the ‘hypercharge chemical potential’ 

[i Y = g\Bo- ( 2 . 6 ) 

It is, like Bq, purely imaginary. 


3 Dimensional reduction 

A useful tool for consistently treating the contributions from the different momentum 
scales at high temperature is dimensional reduction [5J SI 0 IE] • The constant gauge field 
modes (see section [2]) can also be conveniently treated within this framework. Thus the 
computation of the grand canonical partition function is conveniently done as follows: 
In a first step one integrates out hard field modes with momenta of order T. This 
includes all fermion fields because in the imaginary time formalism their (Matsubara-) 
frequencies cannot vanish and are always of order T. The result is an effective action 
containing hi hard; which aside from the zero modes is field independent, and an effective 
Lagrangian Jf so ft f° r a 3-dimensional field theory, and momenta of order gT or less. In 
a second step one integrates over soft modes which are the zero frequency modes with 
spatial momenta of order gT. This yields fi go f t plus an effective Lagrangian for the 
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ultrasoft (p <C gT) fields J?f u itmsoft- When the Higgs mass in 2z? so f t is small compared to 
gT, there are also important contributions from an ultrasoft spatial momentum scale 
smaller than gT , as will be discussed below. After these steps one obtains and from 
that using (I2.3[) . In this way we obtain the grand canonical potential as a sum of 
three parts, 

H = H hard T fi so ft T ^ultrasoft- (3-1) 

In principle it would be possible to treat the constant mode of the gauge fields as 
part of the 3-dimensional gauge held, without introducing the notion of a gauge charge 
chemical potential. Then the distinction between constant and non-constant gauge 
fields would have to be made only when integrating out the soft fields. In such an 
approach the mass term for B 0 would not only contain the Debye mass for the soft 
held, but also a linear and a quadratic term in the constant mode. This point of 
view was taken in [3], For a next-to-leading order calculation it is more convenient to 
distinguish the two as in [2], because the masses for the non-constant modes are only 
needed at order g 2 T 2 , while g 2 T‘ 2 p 1 Lp ~ g A T 2 Bl. Furthermore, in this way we can easily 
read oh the fermionic contributions to fl from [2j. 

4 Hard contributions 

We compute Dhard hi the Standard Model in 4 dimensions. We need the terms of the 
Lagrangian which contain the Higgs held ip, 

= — pd D 2 p — — A ((pV) 2 

- [( h e) a b l a ,LWb,R + (K) ab Qa,L<PUb,R + (Mab Qa,L(fdb,R + h.C.] . (4.1) 

We treat all particles as massless and perform a perturbative expansion in the param¬ 
eters mg, A, hi, and g i: where g 2 and g 3 are the weak SU(2) and color SU(3) gauge 
couplings, respectively. We treat all couplings as being of order g, and ml g 2 T 2 . We 
use dimensional regularization by working in d = 3 — 2e spatial dimensions. Then in¬ 
frared divergences coming from massless propagators vanish automatically. The Higgs 
chemical potential (see (12.61) 1 introduces the following additional terms: 

= fj, v (d T p) - (d T p ] ) ip\ + /4<pV + ZgiTipBoV^T + 2g 2 p ip p ] A 0 p. (4.2) 

Even though we only need an expansion up to order /r^, we find it convenient to include 
the quadratic term in (j4.2j) in the Higgs propagator and later expand the loop integrals. 
Note that there are also /^-dependent vertices whose effects cannot be covered by a 
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frequency shift in the propagator. We will see that the diagrams containing these 
vertices vanish at order /i 2 because the sum integral (14.5j) is zero. 

In the calculation for the hard contributions the following 1-loop sum-integrals 
t p = T E P0 Jp with f p = (2vr)- rf J d d p appear: 

JoM = i H~p 2 ) = - 77 =- ~ 4y + °(a4)> ( 4 - 3 ) 

V 

= = ( 4 - 4 ) 

V 

Here and below we denote p 2 = pi — p 2 , and po = in2iiT + with summation over all 
integer n. The only 2-loop sum-integral which cannot be reduced to products of 1-loop 
integrals is only needed at zero chemical potential, where it vanishes exactly, 


Ji = 


p 2 q 2 (p + q) 2 


P,Q 


= 0. 


(4.5) 


/i =0 


This result has been found to order 0(e) in [9, TO], and to all orders in m- 

We then obtain the following contributions to —fl y, ar r i/H: The leading order is given 
by the 1-loop diagram 


i — —2Jo(pp) — 2 


Vr 2 T 4 2 T 2 4, 


There is also one 1-loop diagram with a Higgs mass insertion 


/T 2 


/+ 


= -2.A!'/«») = —2m 0 -=-s= + 0(/j 


V12 8ir 2 


(4.6) 


(4.7) 


At 2 loops we have the Higgs self interaction, 

2 V A 


} = - 6XJ((pp) 


T 2 ( T 




4 


12 4tt 2 


+ 0(p t 


(4.8) 


The results for the individual diagrams are in Feynman gauge, and we have checked that 
their sum is gauge fixing independent. The gauge fields carry zero chemical potentials, 
and we denote their momenta by q. Their interaction with the Higgs field gives 
1 /-. __d+l 



2 

d + 1 


{g 2 + 3^2) 'A(av)'Mo) 


(. 91 + 3^) 


T 2 

12 


T 2 

12 


87T 2 


+ 0(4) 


(4.9) 
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*- j — jiffl + 3 ^ 2 ) 


2 \ 


Finally, the diagram 

1 


(2 p + qf 


4 W1 ' ^ p 2 q 2 {p + q) 2 

p,g 

= 7(01 + 3^2) [471(^)71(0) — Ji ( n < p ) + 0(/i^)] 

+ 0(/4 


1 

4 l 

^ 2 ■ T 2 f _ 2 


= 4^ + 3^) 12 


12 


iV , ^..4- 

87T 2 


x t www f* — — 7A i ip(fi , i + 3^)72 + O(p) — O(p). 


(4.10) 


(4.11) 


contains the 3-vertices in (j4.2jl which are proportional to p^. Thus at second order 
in p^ we can evaluate the sum-integral with zero chemical potential in which case it 
vanishes, see (14. 5p . The 2-loop contributions above contain symmetry factors 1/2 which 
we have displayed as explicit prefactors of the diagrams. 

All terms of the contributions to fl computed in [2] containing fermionic chemical 
potentials or Yukawa couplings are hard@ Therefore by combining the hard purely 
bosonic contributions computed above with the ones containing fermions from j2j we 
obtain the complete hard contribution as 


12 


VT 2 L 


fl — Q(p = 0) 


hard 


= 6 


1 - 


8vr 2 l 36 4 3 


8tt 2 V 9 3 


tr(/x 


( H 


1 - 


87T 2 

JLfi. 

8vr 2 V 4 
3 


g\ 4^ 

9 3 


+ 


3^2 


tr(^d) 

tr(/i?) 


87T 2 


9 1 


tr(A 


1 + 


mo 

Ji2 


3 fl, , + 3^ 2 

4^\2 x+ ~r~ 

te(K h l)gl - ^ tr ( h l h ug 2 q + KKgl) 


°This is easy to see since the integrals for diagrams with fermions can be written as products of 
1-loop integrals. 
































47T 2 

1 


tr(V4)/4 - + MW) 


47r 2 tr (^^)^ - g^2 tr (/l+/l e /if + /l e /l^e 


<V 


(4.12) 


Here the hi, are the Yukawa coupling matrices (see (14.11) ). The chemical potential 
matrices are determined by the zero mode Bq, or hypercharge chemical potential, and 
by the chemical potentials in (11.11) . 


ha VapY T ^ ^/ I ). ( t • (4.13) 

i 

The T i o are the generators of the symmetry transformation corresponding to the charge 
Qi, acting on fermion type a with a G {q, u, d, £, e}. For example, the generator 
matrices oiB—L are proportional to the unit matrix, with Tb-l ,<j = Tb~l,u = Tb~l,<i = 
1/3 and Tb-l/ — Ts-L,e — ~ 1- 


5 The dimensionally reduced theory 

Aside from the hard contribution f\ard the hard modes also determine the effective 
Lagrangian for the bosonic modes with zero Matsubara frequency, and with soft or 
ultrasoft momenta. The derivation of an effective three-dimensional theory of the 
Standard Model has been done in [8j at zero fi. At order g 2 we need the following 
/^-independent terms@ 

— =^soft, nip=o = ~^Fi jFij + — WijWij 

+ (p Jf 'D 2 (p + + A 3 (v?V) 2 

- l(SiBo) 2 - - i(AA„) 2 - im^Ur (A 2 ) 

- hiy ] (pBl - h 2 </?VTr (Aq) . (5.1) 

For the hnite density effects we also need to include 

- SJ. z? soft = -/i^V - Pi<p*B 0 <p - p 2 ^A 0 ip. (5.2) 

The quadratic scalar operators can be combined, yielding a p v dependent mass j6] H] 
m 3,^ = + m 3 

= m l - hi + T 2 + ~i6^ 2 + 16 ^ + 4 ^ 2 ) ’ ( 5 ‘ 3 ) 

6 The term ip'AoBoip term does not contribute at 0(g 2 ). 
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where h t is the (real) top Yukawa coupling. As discussed at the end of section [3] the 
Debye masses for A 0 , B 0 are only needed at order g 2 T 2 [8], 


m 


D, 1 


m 


D, 2 


N s 5nr\ 9m9 
6 9 ) yi ’ 

2 , N s 5n f \ 2 2 

3 + T + T 92r - 


(5.4) 

(5.5) 


where iV s = 1 is the number of Higgs doublets and nf = 3 is the number of families. 
The couplings are only needed at tree level, 

9i, 3 = 9iT (i = 1,2,3), A 3 = XT, hi = g\y\T, h 2 = J g\T (5.6) 

and also the new parameters in SJf soft , 


Pi ^pipl/ipdl) P2 ^Pip92- 


(5.7) 


In our calculation for the soft contributions we encounter the standard 1-loop inte¬ 
grals 

,3 


T( \ /', / ,2 , 25 2m d T(1 

lo(rn) = / in (A; + m ) — 


him) = 


(k 2 + m?) 


= m 


ll nf 5 

d (4n) d / 2 67t ’ 

«r(i- 


(4:7l) d / 2 ~ - _ 4 ^ + 0(£) 


In the case m = we expand in powers of /i^,, 

/o (m 3 ^) = -^ + ^ + 0(4), 

m 3 ■ tfp , ^7..4' 


The only 2-loop integral we need is PE] 
I(m a ,m b ,m c ) = 


+ 0(/4). 


,.k 2 (k? + m 2 a ) (k| + mg) [(ki + k 2 ) 2 + m 2 ] 


167T 2 


4e 


In 
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m a + m b + m c 


+ 


+ 0(e). 


(5.8) 

(5.9) 

(5.10) 

(5.11) 


(5.12) 


where /i is the MS scale parameter. In the special case m a = m 3i/1 , m b = m e 
{ 0 ,mgi,mD : 2 } and m c = m 3)At it is useful to expand in g 2 


h m 3,^,m,m 3 ^) 


167T 2 

1 

167T 2 


— + In 

4e 


P 


2m 


4e 


In 


3 l/V + m 

P 


P V 


1 

2m 3 + m J 2 m 3 (2m 3 + m) 


(5.13) 

O(pi). 
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6 Soft contributions for soft Higgs mass 


In this section we consider temperatures high enough so that is of order ( gT ) 2 
and positive. At lower temperatures, close to the electroweak crossover, the thermal 
mass squared can be almost canceled by the negative zero temperature ttiq, making m 2 
smaller than 0(g 2 T 2 ). This case will be discussed in section [3 
At 1 loop we have 


) = - 2T/ 0 (m 3iW ) = 2 T 




67T 


47T 


+o«: 


( 6 . 1 ) 


At 2 loops the Higgs self-interaction gives 

2 \ A 


) = - 6AT 2 [h(m 3 ^)\ 
3A T 2 


87 r 2 


(™! -/ 4 ) + 0 ( 4 ). 


( 6 . 2 ) 


Note that the / 1 2 -term has the same parametric form as the one in (14.81) . The sum of 
(16.21) and (14.81) yields the 0(A) correction, that has been computed in [2] by a Higgs 
mass resummation. The interaction between Higgs and the gauge fields gives 

I" (2kx + k 2 ) 2 

/k llk2 ( k i + ™i l/v ) k i[( k i + k 2 ) 2 + m 2 1 


1 /.\ t 2 9 

_ p- j = — (g\ + 3gl) 


3,/i^J 


T 2 r 2 1 

= - -j-(ffi + ^ 2 ) { \h(m 3 ^)] + 4m|^/(m 3) ^, 0, m 3i/v ) j 

V 2 J 2 


32vr 2 


(g'i + 3 gj) 


1 1 

-1-h 2 In 

2s 2 


2m?, 


+ ■■■ 


(6.3) 


1 /' 


2 \ 


1 / 
2 \ 


= ~ tfpT 2 [A(™3,^ , mD,l, "3,^) + 3 A( m s,^ , m D,2, 


g%T 2 

32tt 2 


2 2 


— + 1 + 2 In 
2s 




+3^2 


2s 


+ 1 + 2 In 



2m 3 + m 0i i 

—-— ) 

2m 3 + m D)2 J _ 

3 ,2rp2 ; 


= / i( m 3,^) / i( m ^,i) - 2^ r J i( m 3,^)A(m D , 2 ) 


T 2 


32tt 2 
^T 2 1 
327T 2 2m 3 


(^i' m 3,^m Bi i + 3 j 2 2 ffl 3 ^mD, 2 ) 


+ ZglrriDz) H- 


(6.4) 


(6.5) 
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where we omitted terms of orders other than /^. Adding up all contributions we obtain 
the finite result 


12 


VT 2 


fi(/i) - fi(0) 


= 2/i 


soft 

3m 3 9A 


+ 


+ 


7 iT 47t 2 327t 2 


+ 3.92 C’a] \ + 0(/i, 


4 ) 


( 6 . 6 ) 


with 


Ci 


m D>i 

m 3 


— 1 — 4 In 


( 2m 3 A 

\2m 3 + m Dii y ' 


(6.7) 


After integrating out the soft fields we are left with an effective theory for the ultrasoft 
ones. For soft m 3 the ultrasoft sector contains only the spatial gauge fields. At the 
order we are considering the effective Lagrangian is independent of (i^, so that this 
sector does not contribute to the susceptibilities, and O u itrasoft = 0. 


7 Ultrasoft Higgs mass 

When m\ in (15.ID becomes small, the perturbative expansion used in section |6] can 
break down, which can be seen in (16.5j) where m 3 appears in the denominator. This 
term is of the same order as the soft 1-loop Higgs contribution if \m\\ < g 2 Trri£, g 3 T 2 . 
For such small m 3 it is necessary to include the Higgs held in an effective theory for 
momenta <C gT , which is obtained by integrating out the temporal components of the 
gauge fields. 

First consider H so ft- Since m 3 -C gT we have to put m 3 = 0 in the diagrams in 
section IH1 Then the only non-vanishing contribution comes from the diagram (16.4j) with 
m 3 — y 0. The other diagrams in section [H] vanish in dimensional regularization. Then 
r2 S oft contains an infrared divergence which will cancel against an ultraviolet divergence 
in fiuitrasoft, leaving an order g 2 ln(l///)T 2 /i 2 contribution to H. 

The effective Lagrangian for the ultrasoft fields now reads 

11 2 
-^ultrasoft = - -FijFij + -WijWij - <4D V + 7n 3 ,^V + A 3 (<f^Y (7.1) 

with the parameters [8] 

ml = ml — 4- (3h 2 m D)2 + y v him D y (7.2) 

4tt 

A 3 = A 3 . (7-3) 
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The negative 0(g 3 T 2 ) contribution to results from integrating out the temporal 
components of the gauge fields. It leads to interesting effects depending on how soft 
m 3 is. 

Here we have to distinguish several cases. Consider first m\ ~ g 3 T 2 and positive. 
Then we are still in the symmetric phase. The loop expansion parameter is now g 1 ^ 2 . 
The next-to-leading order (NLO) starts only at 0(g 3 / 2 ) coming from the 1-loop diagram 
(16.ip . and the 2-loop diagrams (16.2p and (16. 3 p contribute at order g 2 . Combining this 
with the soft contribution we find 


12 

VT 1 


fi(/i) - fi(0) 


soft+ultrasoft 



3777.3 9A 3 

7 tT 47T 2 ^ 327T 2 


i + 3g 2 2 C 2 ) 


(7.4) 


with 


Ci 


—1—4 In 



(7.5) 


Note that in this expression we have parametrically ln(7?7£>j/?773) hi(l /g). 

There is another way to obtain (17.41) . Since we are only interested in the O(g^) 
terms we can expand the path integral 


exp( flultrasoft/^"') 


^$ultrasoft exp 


rflrp C/> 

Lt -c«^ u itrasoft 


(7.6) 


to second order in In (17.ip g v only appears in the effective Higgs mass so that 


H(/i) -0(0) 


ultrasoft 


-vt^{^ v ) + 0(4). 


(7.7) 


The expectation value of <p'(p has been extracted from the 2-loop effective potential 

mum, 


(^V) 2-loop 


777 3 T T 2 

2,71 167T 2 


6 A + ( g\ + 3^2) 



(7.8) 


which again leads to (17.4p . 

However, (I7.7j) is also valid when m 3 becomes as small as the magnetic screening 
scale g 2 T of the electroweak theory. I 11 this case the only momentum scale left is g 2 T. 
In a non-abelian gauge theory the physics at this scale is non-perturbative, and the 
loop expansion can no longer be applied, which is the so called Linde problem [ 12] . 
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Nevertheless, the expansion in g (modulo logarithms) still exists, only the numerical 
coefficients in the series cannot be computed by summing diagrams. 

Since the 3-dimensional helds have mass dimension 1/2, and since the only mass 
scale in the ultrasoft theory is g 2 T , we have {(p'tp) ~ g 2 T. Thus the ultrasoft fields 
contribute to hi at order g 2 . A reliable determination of (<p'<p) can only be done by 
lattice simulation of the 3-dimensional gauge plus Higgs system. A recent lattice study 
with mu = (125 — 126) GeV for a SU(2)+Higgs theory can be found in [16]. An 
older but more comprehensive study of the SU(2) theory can be found in [T7] and a 
study including the U(l) gauge fields has been performed in [IS]. Near the electroweak 
crossover (<pV) turned out to be a rather smooth function of the temperature. 

Finally, for negative m 2 the Higgs field develops an expectation value, which in 
presence of chemical potentials for global charges also leads to a non-zero expectation 
value of the temporal component of the SU( 2)-gauge field [4], We have not studied 
this case. 


8 Relation between B and B — L 

To illustrate the use of our results for 11 we compute the relation between the baryon 
number B and B — L in the symmetric phase, which was done in [4] at leading order 
and non-zero Higgs expectation value. First we enforce the saddle point condition 
(12.411 to determine IT as defined in (11.31) . thereby eliminating Bq. Then using (11.41) 
and similarly for (B — L) we express the chemical potentials in terms of B = (B) and 
L = ( L ) which yields a relation 


B = k{B-L). (8.1) 

For m 3 of order gT we obtain using (16.6j) 

_ 4(2 n { + N s ) m 3 24 n f N s 

K _ 22n f + 13 N s + 7 tT (22 n f + 13W) 2 

g\ 236rif - (12 C x - 212)n f W + 75W/ 

+ (22n f + 13W) 2 

g 2 9(12nf — 4(C 2 — l)n { N s + 3N 2 ) 

+ 16vr 2 (22n f + 13W) 2 

g 2 96(8 n 2 + llrifN s + 3 N 2 ) 

~ l&B (22n f + 13W) 2 

h 2 6(6n 2 - 41 rifN 8 - 18A^ 2 ) 

+ lffi^ 2 (22n f + 13W) 2 
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Figure 1: Size of the radiative corrections to k defined in (18. ip relative to the leading 
order result with m# = 126 GeV. The electroweak corrections are rather small, and 
the perturbation series is well behaved. The complete NNLO is dominated by the QCD 
corrections except at the highest temperatures. 
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( 8 . 2 ) 


with the same definitions as in (15.5p and (16.7p . When m 2 ~ g 3 T 2 the result for n can 
be obtained from (18.21) by replacing m 3 by m 3 and Ci by Ci defined in (17.21) and (17.51) . 


The size of the corrections to k are shown in figure [T] over a wide range of temper¬ 
atures. The next-to-leading (NLO) corrections are entirely due to the Higgs, and they 
are quite small. The next-to-next-to-leading order (NNLO) is significantly larger. This 
is caused by the relatively large QCD corrections. When the QCD corrections are left 
out, the remaining NNLO corrections are even smaller than the NLO, indicating that 
the perturbation series is well behaved. We also find that the NNLO Higgs correction 
has about the same size as the electroweak corrections coming from other chemical 
potentials. 

Figure [2] shows a closer look at the most interesting region near the electroweak 
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Figure 2: The ratio of B and B — L at low temperatures with m# = 126 GeV. Shown 
are the LO, NLO and the NNLO result with soft and ultrasoft effective Higgs masses. 


crossover at T ~ 160 GeV. When m 3 is treated as soft, the NNLO corrections diverge 
like l/m 3 when m 3 approaches zero. The perturbation series should be improved at 
small m 3 by assuming m 3 ~ g 3 ^ 2 T and using (17.4[) . It then diverges logarithmically 
when m 3 vanishes. Clearly, the loop expansion breaks down here. However, since (tp'tp) 
is rather smooth when computed non-perturbatively on the lattice, we expect that the 
result for k using (17.71) with the non-perturbative [IS CH IB] should be rather 

smooth as well. 


9 Conclusions 

We have computed the 0(g 2 ) Higgs contribution to the susceptibilities in the symmet¬ 
ric phase of the Standard Model, thus completing the 0(g 2 ) calculation of j2]. Close 
to the electroweak crossover the loop expansion breaks down, and the infrared Higgs 
contributions are determined by the non-perturbative electroweak magnetic screening 
scale g 2 T. Nevertheless, the corrections are parametrically of order g 2 . We have ob¬ 
tained a relation which can be used to determine its coefficient by a lattice simulation 
of the 3-dimensional gauge held plus Higgs theory. We have applied our result to 
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compute the relation of B and B — L. The corrections are small in the regime where 
perturbation theory is valid. Our results indicate that this holds even when perturba¬ 
tion theory breaks down. We find that the QCD corrections dominate except at the 
highest temperatures, and that the corrections are below 5%. 

For leptogenesis our result completes the 0(g 2 ) computation of the washout rate 
|2j. Now two out of three ratesQ entering leptogenesis computations have been obtained 
at this order, the only missing piece being the CP- asymmetry. 
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